Partially locally rotationally symmetric perfect 



FLUID COSMOLOGIES 

Nazeem Mustapha* George F R Ellis| 
Henk van Elst* & Mattias Marklund §1[ 

Cosmology Group, Department of Mathematics and Applied Mathematics, University of Cape Town 
f"^ ' Rondebosch 7701, Cape Town, South Africa 

o ' 

(N| ■ May 16, 1999 

3 



(N 
> 



X 



Abstract 



We show that there are no new consistent cosmological perfect fluid solutions when in an open 
neighbourhood U of an event the fluid kinematical variables and the electric and magnetic Weyl 
curvature are all assumed rotationally symmetric about a common spatial axis, specialising the Weyl 
curvature tensor to algebraic Petrov type D. The consistent solutions of this kind are either locally 
rotationally symmetric, or are subcases of the Szekeres dust models. Parts of our results require 
the assumption of a barotropic equation of state. Additionally we demonstrate that local rotational 
£SJ 1 symmetry of perfect fluid cosmologies follows from rotational symmetry of the Riemann curvature 

tensor and of its covariant derivatives only up to second order, thus strengthening a previous result. 

Preprint: uct-cosmology-99/09 
*q 1 PACS number(s): 04.20.-q, 98.80.Hw, 04.20. Jb 

cr 



1 Introduction 



Our long-term aim is to determine all perfect fluid cosmologies which cannot be invariantly defined by 
the existence of a unique eigentetrad for the fluid shear tensor cr. This consideration is of relevance 
to the equivalence problem of spacetimes (5[ |l7| and its application to relativistic cosmology. The fluid 
shear plays a central role in the dynamics of a generic cosmology If it vanishes, the resultant 
consistent models are highly special and well-studied for perfect fluid matter sources ||. In general, 
cr has two distinct eigenvalues, which may serve to invariantly classify cosmologies on the basis of the 
related eigentetrad. When these eigenvalues coincide (degenerate), but with cr ^ 0, the cr-eigentetrad is 
no longer uniquely defined and the fluid shear tensor thus picks out a distinct spatial direction. It is the 
determination of all cosmological perfect fluid solutions exhibiting this feature that interests us in the 
long run. 

Key symmetries of spacetimes are their continuous isotropies, and cosmologies are either isotropic (and 
then have a Robertson-Walker metric), locally rotationally symmetric ('LRS') (and then are all known 
up to the form of their metric), or are anisotropic (see, e.g., and ]l0| for a discussion of these cases). 
In the case of LRS cosmologies, there is at each event (relative to the family of fundamental observers) 
precisely one preferred spatial direction, and all physical properties and observations are invariant under 
rotation about this direction. It follows that these spacetimes are invariant under multiply transitive 
groups of isometries |2l], |l3j . The question that is interesting from both the physical point of view, 
and in terms of determining the equivalence of cosmologies, is how weak we can make the assumptions 
of rotational symmetry and still determine explicitly the family of cosmological spacetimes involved; in 
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other words, how few physical and geometrical quantities we can make rotationally symmetric, where 
rotationally symmetric means the quantity concerned is either isotropic, or invariant under arbitrary 
rotations about a preferred spatial axis. 

Central to the equivalence problem formalism are the components of the spacetime Riemann curva- 
ture tensor R and its covariant derivatives V . . . VR in a standard tetrad. In Q it was shown that a 
spacetime will be LRS if all tensors algebraically defined by the spacetime Ricci curvature tensor and 
their covariant derivatives up to third order are rotationally symmetric (note that the fluid 4-velocity 
field u is algebraically determined by the Ricci curvature tensor, through the Einstein field equations). 
The ultimate aim is to weaken this assumption by considering perfect fluid cosmologies in which only 
the fluid shear tensor a (and not its covariant derivatives) has this symmetry. All cosmologies that do 
not satisfy this restriction can be invariantly defined through tensor components relative to the unique 
<x-eigentetrad. In the present work we set off to consider a subcase of this more general project in 
the context of an approach that is applicable to investigating systematically all cases with rotationally 
symmetric cr. In detail, we consider perfect fluid cosmologies in which all fluid kinematical and Weyl 
curvature variables are rotationally symmetric about the same spatial axis (the Ricci curvature tensor 
components are automatically so, because of the perfect fluid assumption), thus restricting the cosmolo- 
gies to algebraic Petrov type D and simpler cases. We make no similar assumption about the covariant 
derivatives of these variables. We find all perfect fluid cosmologies satisfying this restriction. 



1.1 Assumptions 

We describe a spacetime by a pseudo-Riemannian manifold M. with a rank two symmetric metric tensor 
field g. It is convenient to decompose the Riemann curvature tensor R into the completely tracefree 
Weyl conformal curvature tensor C a bcd, the Ricci curvature tensor R a b := R c ac b — R(ab) an d the Ricci 
curvature scalar R := g ab R a b according to[] 

Rabcd — Cabcd + R a [c9d]b - Rb[c 9d]a ~ \ -Rffa[c3d]b • (1) 

We will mean by a cosmology a spacetime which satisfies all of the following requirements ^, |l^| : 

• It is a self-consistent solution of the Einstein field equations ('EFE'), relating matter source fields 
which are represented by an energy-momentum-stress tensor T ab to the Ricci curvature tensor R ao 
and its trace R asQ 

Rab - \Rg a b= T ab ■ (2) 

• It is filled with matter energy of some sort which we can represent as a perfect fluid with fundamental 
unit 4-velocity field u (u a u a — — 1). The energy-momentum-stress tensor T ao is then given by 

Tab = fiu a u b + ph ab = T(ab) , (3) 

where h a b ■— g a b + u a Ub denotes the tensor projecting into the rest 3-spaces of observers moving 
with 4-velocity u, and fi :— T a bU a u b and p := ^T a bh ab , respectively, are the total energy density 
and isotropic pressure measured in these rest 3-spaces. 

• Additionally we assume 

(p+p) > , 

so that by (||), together with u uniquely defines the timelike eigenvector of R a b- 

• Observations have shown that the Universe is expanding. We will take this to mean that the 
(isotropic) expansion of u, described by the fluid expansion scalar O, is positive: 

6 > . 



1 When referring to tetrad components in the spacetime, we use Latin indices from the first half of the alphabet: a, b, c € 
{ 0, 1, 2, 3 }; in the spatial sections we use Greek letters: a, f3, 7 £ { 1, 2, 3 }. For the local coordinate spacetime description 
we use Latin indices from the second half of the alphabet: i, j, k 6 { 0, 1, 2, 3 }. Round brackets denote symmetrised indices 
and square brackets denote skew-symmetrised indices. 

2 Throughout this work we will employ geometrised units characterised by c = 1 = 8irG/c 2 , and we set the cosmological 
constant A equal to zero (the latter implies no loss of generality, as A can be effectively included by suitably redefining T a t,). 
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The Weyl conformal curvature tensor may be decomposed relative to u into its symmetric tracefree 
'electric' and 'magnetic' tensor parts, E and H, respectively, according to || |8| 

E ab ■= C cedf h c a u e h d b u f = E {ab) (4) 

H ab := (-hVceghC 9h df )h C aU e h d b uf = H {ab) , (5) 

with the completely skew spacetime permutation tensor r\ a bcd specified by r\ abc d = V[abcd] an d 770123 = 1, 
77 0123 = — 1. The second Bianchi identities differentially relate components of the Riemann curvature 
tensor: 

V [a i? bc]de = . (6) 

As well as entailing the matter equations of motion VbT ab = = V b (R ab — \R g ab ), this relation provides 
evolution and constraint equations for E and H when employing the 1 + 3 decompositions (0) and (ra) 




1.2 The problem: symmetry assumptions 

The question we answer here is: 'Under what conditions are all of the fluid acceleration, vorticity and 
shear and the spacetime electric and magnetic Weyl curvatures either rotationally symmetric about the 
same spatial axis or isotropic, but the spacetime itself is not?'. A dynamical tensor field is rotationally 
symmetric if there is a degeneracy in that tensor quantity in terms of its eigenvalues, but it is not 
isotropic (not all eigenvalues are zero). The justification for this terminology is based on the fact that 
a spacetime is LRS if all tensor quantities, as well as their covariant derivatives, are either isotropic or 
rotationally symmetric about the same spatial axis, with at least one not being isotropic J?], Our 
ultimate aim is to be able to uniquely classify all perfect fluid cosmologies which have degenerate fluid 
shear eigenvalues. To this end, we employ the methods developed in ]7j which introduce as a reference 
framework an orthonormal tetrad ('ONT') field { e a } where the timelike member e is identified with the 
fluid 4- velocity field u ('1 + 3 decomposition'). For the spatial rotation coefficients we use the notation 
suggested in || [18], [ll| and utilised to great effect by Wainwright and collaborators (see p5j for a survey). 
We have thus written out for a perfect fluid matter source the evolution and constraint equations 3 for 
the tetrad commutation functions 7 as well as for E and H, which are obtained from the Jacobi, Ricci 
and second Bianchi identities, and specialised to a rotationally symmetric fluid shear tensor er. These 
are given in full in appendix |A|. 



2 PLRS perfect fluid cosmologies 

The effect of restricting the geometry of a spacetime is generally expressed by setting certain geometri- 
cally defined tensor components to zero. The evolution equations for these quantities now become new 
constraint equations. These constraint equations must be preserved along the fluid flow lines u to be 
consistent. In general this does not happen without further constraints on the dynamics. Of course, 
the test for preservation then has to be repeated for these new constraints. We list all perfect fluid 
cosmologies that provide solutions to the EFE which have the tensor fields er, E and H as well as the 
vector fields u and u), but not their covariant derivatives, all rotationally symmetric about the same 
spatial axis. In our investigation we find that all consistent cosmological solutions to the EFE satisfying 
this criterion, bar one — the Szekeres dust solutions — , are in fact LRS spacetimes in the definition of 
Ellis and Stewart and Ellis pl| . We shall generally take it / 0, as for cr — we will see that the 
presently considered types of cosmologies typically specialise to the well known spatially homogeneous 
and isotropic Friedmann-Lemaitre-Robertson- Walker ('FLRW') cases. 

Definition 2.1 We assume that in an open neighbourhood U of an event on an expanding perfect fluid 
spacetime manifold M. 

((722 — 0-33) = 012 = 023 = 031 = , 

(£"22 — £33) = Eyi — E23 = E31 = , 

3 When eo = u has non-vanishing vorticity this terminology is doubtful because then the local rest 3-spaces of neigh- 
bouring fundamental observers do not mesh to form everywhere spacelike 3-surfaces. 
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(H 2 2 — #33) = H12 = H23 = H31 = , (7) 

ii2 = U3 = , 
oj 2 = W3 = . 

That is to say, we have set all of cr, E and H simultaneously rotationally symmetric about the same 
spatial axis. In addition, the vectors u and uj are aligned with this axis of rotational symmetry. We call 
all cosmological spacetimes in which there exists a tetrad { e a } with respect to which these restrictions 
are satisfied partially locally rotationally symmetric ('PLRS') cosmologies. Those PLRS cosmologies that 
are not LRS will be referred to as strictly PLRS. 

Remarks: 

(i) Using a er-eigentetrad { e a } here is equivalent to using an E- or H-eigentetrad. 

(ii) If either of the vector quantities u or uj do not point in the invariant spatial direction defined by 
the degenerate tensor quantities, then they define an additional invariant spatial direction. The tetrad 
{ e a } may then be invariantly defined by, e.g., identifying another of its members with this new direction. 
And thus, for the situations we want to start from to be PLRS, we must have that ii2 = «3 = 0; that is, 
u || ei. Similarly, we must have w 2 — ^3 — 0, corresponding to uj \\ e±. If these conditions did not hold, 
this would exclude for non-zero u and uj the possibility of the spacetime being LRS and thus compromise 
the presently proposed notion of partial symmetry. 

(iii) The present setup arbitrarily adapts to the spatial ei-axis. However, this is only a matter of 
convention and by a cyclic permutation of indices 1 — > 2 — > 3 — > 1 one can easily adapt to any of the 
other spatial axes as well. 

(iv) Our assumptions reduce the perfect fluid cosmologies we consider in this article to algebraic 
Petrov type D; a classification scheme for solutions with this property based on the Newman-Penrose 
formalism has been suggested by Wainwright in (24). 



With Definition [2.1| in place, the momentum conservation equations for all PLRS cosmologies reduce 



i n U to ( P7| ) - (|99|). Moreover, combining the H-constraint equations (108) with (IPS) and ( 110 ) with 



(111), we find that all PLRS cosmologies require in U that 

(n 2 2 - n 3 3) (ofi +|wf) = n 23 (o^ + § uj\) = 0; 
which means, as we assume that eruj ^ 0, thatQ 

("22 - "33) = "-23 = . (8) 



The related evolution equations obtained from (|77|) and (|78|) as well as (70) and (fnl), respectively, now 



reduce in IA to constraints on the e.2- and e3-gradients of the Fermi- rotation variables Q2 and O3; namely 

(e 2 + 2n3i)(02)-(e3-2n 12 )(n 3 ) = (9) 
(e 2 + 2n 31 )(r! 3 ) + (e3-2n 12 )(fi 2 ) = 0. (10) 

We proceed to check the consistency of the cosmological PLRS subcase of the EFE by mainly com- 
puting the time evolution of all new constraints. In the process of doing this we will fix the remaining 
tetrad freedom and thus invariantly classify the solutions. For each of the cases we consider below we 
will check the transformation behaviour of the commutation functions and other tensor quantities and 
use this to fix the freedom conveniently. 

We now turn to the issue of choice of spatial triad { e a }. As e! is presently a uniquely defined vector, 
ei is fixed. This will thus have an invariantly defined direction, say X. For any given tetrad choice, the 
components of this fixed direction in the spatial e 2 - and e 3 ~directions are X ■ e 2 and X • e 3 , respectively. 
We have the freedom to set one of these components (or any other quantity which does not behave like a 
scalar under a rotation) to zero in U by rotating the spatial triad { e a } in the e 2 / e 3 -plane. Alternatively, 

4 This restriction on n a p suggests that for PLRS cosmologies according to Definition 2.1 no consistent solutions to the 
EFE exist that contain gravitational radiation; it is the transverse (with respect to the spatial ei -direction) components 
(«22 —^33), fi23 r _(f r 22 — C33) and CT23 which typically form the connection characteristic eigenfields associated with the local 
null cones (cf. |14|). If we relax the PLRS restrictions so that E is not rotationally symmetric, we find that this feature 
still holds. 
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as can be seen from the transformation property of we may choose to set fii = in U which fixes 
the eo-gradient of the rotation angle tp; and then for example we choose either 77 22 or 71,33 m the 3-space 
by fixing the ei-gradicnt of ip. 
For ease of notation we will use 

a :— i + a aa (no summation) . 



2.1 General perfect fluid 

This is the most general case to be considered here where we will take throughout that li^O and u> ^ 0. 
For u> = see section 2.2, for u = refer to section 2.3. We assume a perfect fluid, but leave the equation 
of state unspecified. 



2.1.1 Constraint analysis 

We proceed by fully fixing the tetrad freedom by choosing e3 orthogonal to the projection of the fixed 
vector X in the e 2 / e3-plane. Hence, we rotate the spatial triad such that 

fl 2 = e 3 • ei = -ei • e 3 = . (11) 

Under the assumptions of PLRS symmetry (0), which imply (||), and with the above tetrad choice, 
the E-equations (Q) and ( p5| ) yield the new constraints 

e 2 {H n ) = (12) 
(a 2 -n 31 )H 11 = 0, (13) 

while from the H-equations (^) and the new constraints 

e 2 (/i) = 3(a 2 -n 31 )E 11 (14) 
e 2 (E n ) = ±e 2 (/z) (15) 

arise. The algebraic condition (jl^) suggests that we distinguish between two subcases according to 

A] (o 2 - n a i) = or B] H n = . 



A] (a 2 — 77,31) = 0: So now (a 2 + n 3 i) — 2 a 2 . Fr om the cj-equation (104) we see that e 2 (ui) = and 
from ( |l4|) we have e 2 (/j,) = 0. The <r-equation ( |l02| ) then shows that also o\\ — must hold, providing 
a split into further subcases according to 



Al] ^3 = or A2] cr n = 



Al] SI3 = 0: This has the implications from ( |l03| ) that es(ui) = 0, which implies from (105) that 
(d3 + 77,12) = 0, as we assumed ii\ ^ 0. We note that the following e 2 - and e3~gradicnts of certain 
quantities must vanish: from ( |l2] ) and (p7|) we get that e 2 (Hn) = e3(i?ii) = and from j93| ) - ( ^6|) we 
find that e 2 {E\\) — ea(Eii) = and e 2 {^i) = e^i) = 0. Then (f72|) and (ff^) reduce to 



e 2 (6i) = e 3 (6i) = 



(16) 



Now we check the propagation property of the vanishing e 2 ~ and e3~gradients of the energy density fi. 
We use the commutator relations ([39]) and ( J60| ) operating on fj, and the energy conservation equation ( |64| ) 
to show that e 2 (/j) = e 3 (/j) = if (/i +/j)e2(6) = (fi +p)e3(9) = 0, since the momentum conservation 
equations (^8|) and (|9^) must hold. Thus, if we want to stick to purely cosmological solutions as we defined 
them in the introduction, we must have e 2 (0) = 63(G)) = 0. Now this means that e 2 (tt>i) = 63(^1) = 0, 



which we obtain from the Ricci identities ( p_13[ ) and (115), suitably combined with (|16|). This result 
is crucial because now we can show that the solutions contained in the present PLRS subclass are not 
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cosmological ones. To do so we find the commutator (61) most useful. Wc first apply this relation to W\, 
yielding 

2 6 2 wi - cti nil = , (17) 



where we have used the w-equation (p9|), and the constraint on ei(wi) given by (106). We then apply 
the commutator (plj) to the energy density //, using (p3), and find that 



2 ui (fx + p) — nu ei (fj,) = 



(18) 



Finally we apply the commutator (|6l| ) to the electric Weyl curvature component En and substitute from 
©, (H|) and © to get 



2 2 wi [{n+p) -3 E n ] +3ai n u J5 n = 



(19) 



We now use (|l7|) in (|ll]) to get 62 u>\ (fi + p) = ©2 = 0, as we assumed (/i + p) > and Wi 7^ 0. 
Substituting back into (|l7| ) now give s ai Tin = 0. We argue that this necessarily means nu = 0. If 
instead we started from a\ — 0, then ( |llO| ) shows that n\\ u>\ = => n\i = 0; hence, n\\ = in any case. 
But then (116) shows 0i = 0, and thus, since already 02 = 0, we find with the Raychaudhuri equation 
(||) that in U 

= 0, (20) 

violating the premise that our cosmology be an expanding one. We conclude that there exist no cosmo- 
logically viable solutions in the present PLRS subclass. 



A2] (T11 = 0: From the &— and w-equations (103) and (105) we must have 

(a 3 +n X2 )ui -CI3UJ1 = , 
and from the second Bianchi identities (|8^ ) and (j87l), 

(a 3 + n i2 ) Hn + f2 3 En = . 



(21) 



(22) 



Combining these two equations, we get that cither (a 3 + 7112) = £! 3 = which would be a subcase of Al] 
and is thus non-cosmological, or 

u-i E n + at H n = . (23) 



An important algebraic relation is given by the H-constraint (107); that is 

2 (iii + ai) lui + Hu = . 



(24) 



The tetrad choice employed has the e ffect of gen eral ly eliminating the gradients in the e2-direction of 
important scalars. In particular, from (114) and (115) we get 



Now, from ( |104| ) we get that 
and from ( |130| ) 

From (p3l), suitably combined with 



e 2 (o>i) = = e 3 (0) . 

e 2 (ui) = , 

eatffn) = . 
we get 

e 2 (£; 1 i) = e 2 (M) = 



If we now take the e2-gradient of equation ( p4| ) and substitute from the equations (|2E 
we get the useful result 

e 2 (ai) = . 



§ and (t 



A relation for e 3 (wi) is provided by (112) and ( |113| ): 

e 3 (wi) + (a 3 + m 2 ) wi = 



(25) 

(26) 
(27) 

(28) 

©. 
(29) 

(30) 
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We get relations involving the e 3 -gradients of both En and \x from combining (|95|) and (|96|) suitably, 
thus yielding 

e 3 (£ai)-ie 3 (/x)- 3ft 3 #11=0 , (31) 

e 3 (p) - 3 (a 3 + n 12 ) E n +3ft 3 #ii =0 . (32) 

The commutators provide vital information here. We find from ( |6l| ) acting on wi that 

- | 9 lui + ai nn = 2 a 2 (a 3 + n i2 ) , (33) 

where we have used the <i>-equation (j69|), the constraints on the gradients of u>\ provided by (EOS), (25) 



and (j3Cl), and then substituted from (|116|) with (117) appropriately. We find from the commutator (51) 



acting on the magnetic Weyl curvature component H\\ that 



3 Hn (-36^1+01 nn) — "-11 (" + P) w i — 12 a 2 -En f2 3 = , 



using a relation for e 3 (iJn) provided by (|87|), noting the constraint ( |27| ) on e 2 (iJn), the constraint on 
ei(iJn) given by ( |129| ), and then using the H equation (|88|). We also needed the first part of ( p8| ) and 
(|l0|). Combining the above with (33) we get the useful result 



nu (p + p) ui = - 18 a 2 ^3 En , (34) 
where we have also used (p2[). We now take the commutator (p3[) operating on cji and substitute from 



(p5j), (26), (pj|), the vorticity constraint equation (106), and ( |30| ) into the resultant expression and we 
find that n 33 (a 3 + n i2 ) Wi = 0. Now if (a 3 + ni 2 ) = 0, then from ( pl| ) we must have fl 3 = 0, and thus 
this would be a class already dealt with in section A2]; those were non-cosmological. So we conclude 
that n 33 = 0. We now take the e 3 -gradient of (^) and get the key result 

uie 3 (/i)=0 =► e 3 (/i)=0 (35) 



by using, in addition to (|23| ), equations (|22|), pQ ) - (|32|), ( J105[ ) and (|87|). We proceed to check the 
consistency of ( |3~l| ) and ([32|) in the combined form e 3 (£n) — 3 (a 3 + ni 2 ) En = 0. We propagate this 
along u and find 

^n 11 [{a 3 +n a )H 11 -Sl 3 E 11 ]=0, (36) 

where we have used the E-equation (|79|),the constraints on e 2 (_Bn) and e 2 (/x) given respectively by ( p8| ) 
and (p5f), the constraints on e 3 (ifn) given by (|S7"|), and the Jacobi identity ( f73| ) together with ( |114[ ). The 
evolution along u of the e 3 -gradient of En is obtained by applying the commutator (pm to En and using 



the equations ( |79| ) with (121). Now we may rewrite (|3(j) by using (p2|), obtaining nn (a 3 + n i2 ) i?n = 0. 
We show that this means that nn = 0. If not, then (a 3 + ni 2 ) Hn = which means that f£ 3 En = 
(from (|2^)), which then contradicts (|34|). 

We have seen here that nn = is required. But now, again from ([34]), we must have a 2 f2 3 Bn = 0. 
So either a 2 or 3 £n vanishes. We consider these possibilities below. 



If Q 3 En = 0, then from (|22|) we require (a 3 + ?ii 2 ) Hn = 0. If now Hn = 0, then ( p2[ ) tells us 
that these solutions are not cosmological since they require in U (u + p) u>\ = ^> (u + p) = 0. 
And if (a 3 + ni 2 ) = 0, then from ( |33| ) we get in IA 9cji = 9 = 0, and we are clearly in the 
non-cosmological realm again. 



If, on the other hand, a 2 = 0, then again from ( p3[ ) we must have in U 9 u>\ — 9 = 0. So none 
of the solutions here are of relevance to us. 



B] Hn = 0: Immediately we see from ( |S2"| ) that (p + p) + 3 En = 0. Also, from (|88|), we get that 
"-11 En = 0. Now if En — 0, then (/i +p) = which is not allowed. So nn = 0. We also get from ( p6| ) 
that fi 3 i?n = 0, and once again we deduce that since En = (p + p) = 0, it follows that we must 



have fl 3 = 0. Moreover e 3 («i) = from (h03); and e 2 (Mi) = from (102). This tells us from (104) that 



(a 2 — n 3 i) = and from (105) that (a 3 +ni 2 ) = 0. And now from (117) we get that 9i u>\ — 9i = 



We get from (|93|) and (|9j) that e 2 (E 11 ) = e 2 (u) = and from (g5|) and (|9J) that e 3 (£ai) = e 3 (u) = 0. 
If we now take the above two relations for the e 2 - and e 3 -gradients of /i and substitute them into the 
commutator ( |6T|) applied to fi, using (|6j), we get & uji (fi + p) = 0, and so in U 

9(n+p) = 0; (37) 

in other words, there are no solutions in this PLRS subclass that are of cosmological interest. 
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2.1.2 Summary 

We conclude that there exist no rotating and accelerating perfect fluid cosmologies which are PLRS 



according to Definition 2.1 



2.2 Irrotational accelerating perfect fluid 

These models have u = 0. We assume that u ^ 0. An immediate implication for th e com mutation 
functions n a p, in addition to (||), is that nn = 0, from (p9). Now from the H-constraint (107) it follows 
that Hn = 0, reducing the H-equation ( jS8[) to a trivial statement. A useful point of departure is provided 
by the E-equations (|84|) and (^(]). That is, En = SI3 En — 0. A brief argument below will show that 
this means that in U 

n 2 = il 3 = . (38) 

• The argument goes as follows: if En = 0, then from ( |79| ) we have (// + p) an = an = 0. 
We find that the following gradients vanish: ei(/i) = = ei(6) from ( [126[ ) and (110), respectively. 
Now from the commutator relation (|5^) acting on /i we get [fx + p) =0; that is, this is a non- 
cosmological subcase. In this little argument, we have employed the assumption that the matter 
fluid has a barotropic equation of state, p = p(/i). 



So we must have (|38| ) holding, which then implies from (102) and (102) that 



e 2 (wi) = 

In turn the effect of the above is that (02 — n.31) 
So we may write (02 + 1131) = 2a 2 and (03 — ni 2 ) 



e 3 (ui) = . (39) 
"12) = which derives from (104) and ( |105| ). 



(03 H 
2a 3 . 



2.2.1 Tetrad choice and constraint analysis 

From the above we see that we are free to employ along u a Fermi-propagated cr -eigentriad { e a } by 
setting fli = eo(y) via a spatial rotation by an angle ip in the e2 / e3~plane so that £l[ = 0. We note the 
following consistency conditions: 



nn 

(o 2 - n 3 i) 
(a-3 + n 12 ) 



(from (76)) 

(from @ and fljlg)) 

(from dH) and (Q) 



We can further use the tetrad freedom on a 3-surface x° — c° to set 7133 = 0. This condition is then 
preserved as we can see from ([78]). We can also set 03 = on a 2-surface x° = c°, x 1 — c 1 . This we can do 
with impunity since firstly eo(a3) = — 82 03 from ([75]) and ( |115| ). Secondly 61(03) = a% 03 from (118) and 
(119). We conclude 03 = everywhere. We may summarise: the only non-zero commutation functions 
are u±, Oi, 02, a\ and a 2 - Of these quantities, only a 2 does not have its e2 and e3 gradients vanishing; 
as can be seen from (|39"|), (112) and (113), (114) and ( |115| ), (119) and (121). We may now proceed to use 
the remaining tetrad freedom to set 63(02) = on the line x u = c°, x 1 = c 1 , x 2 — c 2 . This we can do 
by observing that the eo~, ei- and e2-derivatives of this quantity are conserved respectively by applying 
the commutators @, @ and @ to a 2 . We also need @ and ( pdl| ), and ( |i20| ) and ( |i22] ) to see 
this. Finally we set 62(02) = at an event x % — c 1 . This can be done because firstly, the eo-derivative of 
e 2 (a 2 ) is driven by a multiple of e 2 (a2) (from (|59|) operating on a 2 and (f74|)); secondly, the ei-derivative 
of 62(02) is driven by a multiple of 62(02) (from (p5|) operating o n a 2 and (|l20| )); thirdly, the e2-derivative 
of e 2 (a2) is driven by a multiple of e 2 (o2) (from taking e 2 of (|l22j )); lastly, the e 3 -derivative of e 2 (a2) 
vanishes (from ([H]) operating on 02). The remaining commutation functions have their gradients in the 
e 2 I e3~plane vanishing, and the only equations remaining constrain quantities in the invariant spatial 
ei-direction (from (11C), (123) and ( |l20[ )). The now algebraic relation (122) determines £11 in terms of 
H, say. The curvature variables which remain non-zero in U are p, /1, and En- We find that all of these 
variables have their gradients vanishing in the e 2 / e3-plane. This is obvious from (^8|) and (^9[), and 
( p3| ) combined with ([m]), as well as (|95|) combined with (^). The remaining constraints on the various 
ei-gradients are given by (^7j) and (126). The consistent solutions in U we so obtain are the LRS class 
II solutions of Stewart and Ellis lElJ . 
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2.2.2 Summary 

Having assumed in this subsection a barotropic equation of state, p — p(fx), we conclude that the only 



irrotational accelerating perfect fluid cosmologies that are PLRS according to Definition 2.1 are the 
expanding solutions in LRS class 11 of Stewart and Ellis Q] (see also Jl3j and 

2.3 Rotating dust 

These models have p = 0. From part of the twice-contracted Bianchi identities (momentum conservation 
equations), e a (p) = =>• u a = 0. We here assume that u> ^ 0. Since we are dealing with dust, we note 
that now our tetrad choice is that of Ellis in 0J^| However, because we demand 6 > 0, we will refrain 
from identifying his rotating dust solutions of LRS class I which require O = 0. We immediately note 



from the Jacobi identities (104) and (105) that in U 

n 2 = n 3 = o . 

2.3.1 Tetrad choice and constraint analysis 

We are free to propagate the spatial triad { e Q } along u as anti-rotating by choosing (oji = eo(ip), 

i.e., (uji + fli)' — 0. We can further use the tetrad freedom in U to set n 33 — by choosing ei(ip) such that 
ei((p) = — 7133. This quantity is conserved, as we can see from (f78|). But now we may not as yet proceed 
with a further tetrad specification as in (7j, where (03 — 7112) = on a 2-surface x° — c°, x 1 = c 1 , because 
this would constrain the geometry. In particular, it actually requires e2(o>i) = for (03 — 7742) = to 



hold, using (75) and (115). So we proceed by leaving the tetrad freedom unfixed for now and see what 
the implications are from consistency checks. From setting a to be rotationally symmetric we do not 
get any immediate constraints. But from setting E to be rotationally symmetric we get new constraints. 
Specifically we get 

(ci2 — 7131) Hu = (from the E-equation ( |34"| ) combined with ( |130D ) 
(03 + 77,12) H11 = (from the E-equation (|S6|) combined with (|131|)) . 



So naturally here we have a split into 

A] H n =0 or B] (a 2 - 7731) = (a 3 + 77,12) = . 

A] Hu = 0: For Hu to vanish, we must have from the H-equation ( |S8| ) that Tin En = 0. 

• A brief argument now shows that we are not interested in En = 0. It goes as follows: if En = 0, 
then from d73) we have /ion = =>• en = 0. For consistency we now require from the <r-equations 
( |65| ) and (|66|) that lui = 0; that is, this case is dealt with elsewhere. In fact this leads to FLRW 
solutions as we have already noted. 

Then we conclude that Bu /0 must hold and hence fin = 0. This has the immediate implication 



from the H-constraint equation (107) that a\ — 0. Now we get from the second Bianchi identities the 
following new constraints 

3 (oa - nai) £11 - e 3 (/*) = (from (93) and @) (40) 
3 (03 + 7712) E n - e 3 (/j) = (from © and ©) . (41) 



We propagate (J40j) along u twice, using the necessary evolution equations. We use (|72j) , (|79[), and fl59|) 
operating on /j with Q), to get 

wi (a 3 + ni 2 ) (3 £■!!-//)+ 2 //e 3 (o;i)=Q . (42) 



We now need fl69|), (f73j), ([79[), (64), and (BOJ) operating on uji, to get for consistency of (f!2j) that 



- u)\ (a 2 - 7731) (3 E n - /j) + 4/je 2 (wi) = 

6 Described in his Theorem 3.1. 



(43) 
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We propagate ( fill) along u twice, using the necessary evolution equations. These are ((73|), (|79|), and ( |60|) 
operating on [i with (pj), which thus yields 



u\ (a 2 - n 3 i) (3 .En - //) + 2/ie 2 (wi) = 



We now need (|69|), (|72|), (|79|), (|64| ), and ([59|) operating on L)\ to get for consistency of (|4J) that 



cji (a 3 + m 2 ) (3E?n - m) - 4/xe 3 (wi) = . 
We form linear combinations of the above four constraints to facilitate our task at this point. 



(44) 



(45) 



D) 
3)) 



/ie 3 (wi) = (from§- 

/Lte 2 (wi) = (from@ + 

h+«i2)(3£n^) = (from 2 x © + @) 

(03-n 81 )(3J?ii-M) - (from @ -2x @) . 



e 3 (wi) = 
e 2 (wi) = 



• A brief argument now shows that 3 En — fi = is not applicable. It goes as follows: we note 
that e 2 (tJi) = e^(uji) = 0. So from the commutator ( |6l"|) acting on o>i and incorporating the 
adequation ( |69| ) into this, we get that 9 2 = 0. The consistency of this requires from (^) that 
3 En — /x+ 3wf = 0, and if 3 En — fi = 0, it must necessarily follow that u)\ = 0, which we 
excluded. 

So we must conclude that (a 2 — 72,31) = (a 3 + ni 2 ) = 0. This is severely restrictive. We get from the 
commutator (^lj) acting on uj\ that 2 = 0, using the cj-equation (|6^). And now from ( |116| ) we must 
also have Oi = 0. These last two results, in particular, imply that in IA 



6 = 0, 

which means that spacetimes in this PLRS subclass are non-cosmological. 



(46) 



B] (02 - n.31) = (a 3 + n u ) = 0: 



Firstly we note from (112) and (114) that 
ea(e0 = esCGi) = . 



3) combined with 



(47) 
2), and 



The critical constraints are obtained from the second Bianchi identities: 
( p5| ) combined with ( |96| ) once again. They read, respectively, 

e 2 (Sn) = e 2 (/i) = 
e 3 (Sn) = e 3 (/i) = 0. 

We proceed to check preservation along u of e 2 (/z) = e 3 (/i) = by using evolution equations obtained 
from the commutator relations ( |59| ) and (|60|) when acting on /i. We also require from this the relations 
given by @, (11121) , (|Tll ), ([111 ) and (|U5T We get that 



eo(e 2 (/i)) 
eo(e 3 (/i)) 



(30-11 
(2 011 



|e)e 2 0*)-2e 3 (a; 1 ) 
§ 9) e 3 ((i) + 2e 2 ( Wl ) 



e 3 (wi) = 
e 2 (u>i) = 



Taking these results and putting them back into (112), ( 114 ), (113) and (115), and recalling (J47]), we have 
e 2 (Oi) = e 3 (0i) = e 2 (0 2 ) = e 3 (0 2 ) = 0. We note now from (13C) and (131) that it is apparent that 
the gradients of H\\ are also degenerate in th is fa shion: e 2 (iJn) = e 3 (i?n) = 0. We check preservation 
along u of e 2 (i?ii) = e 3 (£ u ) = by using (126) and the commutators ( p39| ) and (|60| ) when acting on 
E\\. The results are that 



e (e 2 (£ii)) 
eo(e 3 (J5n)) 



§ fin e 2 (nn) 
§ Hu e 3 (nu) 



If now Hn = 0, we have a case already dealt with in A]. So we must therefore conclude that e 2 (?in) = 
e 3( n n) = 0- Now here again crucial algebraic constraints are obtained from the commutator (|6ll). Acting 
on wi we get 

9 2 ^i = ai nn (48) 
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from the adequation ( |69[) and the ei (tji )-constraint fllOq ). Acting on nn we get 

nn (o-nWi +ffliriii) = (49) 



from the evolution equation for nn ( |7q ) and the ei(nn) -constraint (116). We show that we do not get 
any relevant solutions here. We start with equation ([i"9|). 



We first show that nn = leads to trivial solutions. If nn = 0, then from (116) we must have 
0i = and then from (jig) this means that 9 2 — and so this cannot be cosmological because 
consequently in IA 

9 = 0. (50) 



• If, on the other hand, nn ^ 0, and instead we have from ( |49] ) that o~n u>i + ai nn = 0, then we 
immediately get from ( ^8| ) that 9cl>i = 0; that is, in U 

9 = 0. (51) 

So these too are not cosmological solutions. 
2.3.2 Summary 



We conclude that there are no rotating dust cosmologies which are PLRS according to Definition 2.1 



2.4 Irrotational dust 



We shall now consider irrotational dust cosmologies under the PLRS restrictions of Definition 2.1: that 
is to say, cosmologies where uj = and p = => u = 0. We can show that the only consistent irrotational 
dust solutions of the EFE which have er and both E and H rotationally symmetric in the same plane are 
known sol ution s ,p| F or ex ample, the case er = leads to the FLRW dust solutions 0. Proceeding, we 



find from fllOfl ) and ( pL02|) that in U 

n 2 = n 3 = o 



2.4.1 Tetrad choice and constraint analysis 

We note that for the models we are interested in here, we again choose the cr-eigentriad { e a } to be 
Fermi-propagated along u by setting fli = e (ip) via a spatial rotation by an angle ip in the e 2 / e 3 -plane 
so that Sl[ = 0. We may also choose n' 33 = ei(ip) = — 7133, since now the evolution equation ( |78| ) for 
n33 becomes involutive. We will show that the present PLRS subclass recovers the well known Szckcres 
dust solutions From the E- equation ( |84| ) we find e 2 (-Hn) = | (a 2 — n 31 )iJ n , while (p5|) gives 

^(-STli) = 0, leading to the more useful algebraic result 

(02 -n 3 i)Hn =0 . (52) 

Then from ( ^6| ) we find e 3 (7?n) + | (a 3 + n i2 ) i?n = 0, while j87| ) gives e 3 (_ffn) = 0, now yielding the 
algebraic result 

(a 3 + m 2 )Hn = . (53) 

Now the commutator ( |6l| ) is helpful at this point. If we operate on Hn, we find that nn ei(iJn) = 0, 
which subdivides the class into 



A] nn =0 or B] ei(JJn) = . 

7 We may, with relative, ease generalise this result to the situation where only <r and H are rotationally symmetric. 
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A] nn = 0: It is easy to see from ( |76|) that presently rin = is conserved along u. Hence, the H 

constraint (107) gives Hn = 0, which the H equation ( |S8|) preserves; thus H = 0. Now we look at (pif), 
which reads 

e 2 (£ u ) = ±e 2 (/i) , (54) 



and we use (|79|), ( |ll3[ ), (|64|), and the commutator ( |59|) acting on En and /i, to check the time evolution 
property of this last relation: 

e (e 2 (£ u ) - ±e 2 (,i)) = -49 2 [e 2 (# n ) -|e 2 (/i)] . 

This is solved by e 2 (i?ii) = |e 2 (/_i), and we may now use ( |127| ) with the above to get e 2 (En) — 
| (a 2 — 7131) E 1 !! — i e 2 (/j), which, in this manner, shows the consistency with (p3|). Now we look at (p6|), 
which reads 

e 3 (£n) = ±e 3 ( M ) , (55) 

and we use (|79|), ( |115[ ), (|64|), and the commutator ( |60| ) acting on En and /j, to check the time evolution 
property of this last relation: 

e (e 3 (£ii) - |e 3 0*)) = - 4 6 2 [e 3 (£ u ) - ie s (A*)] • 



This is solved by e3(_En) = -|e3(/i), and we may now use (|l28| ) with the above to get e3(_En) = 
| (as + ni2 ) En — \ which, in this manner, shows the consistency with (|9^). We have from 

(121) and (HE) that along u, and, by substitution into the relevant commutators, that everywhere 
e 2 (ai) = e^(ai) = 0. We may at this point use some of the remaining freedom in the commutators to 
(on a 2-surface x° = c , x 1 = c 1 ) set e3(a 2 + n^i) = e 2 (a3 — ni 2 ). This can be done because firstly 

eo(e 3 (a 2 + n 3 i) - e 2 (a 3 - n i2 )) = - 26 2 [e 3 (a 2 + n 3 i) - e 2 (a 3 - n i2 )] , 

using the commutator ( |60| ) on (a 2 + n^i) and the commutator ( |5Sj ) on (a3 — ni 2 ). Secondly, 

ei(e 3 (a 2 + n 3 i) - e 2 (a 3 - ni 2 )) = 2a x [e 3 (a 2 + 7i 3i ) - e 2 (a 3 - ni 2 )] , 



using the commutator 
Thus 



on (a 2 + 7i3i), the commutator (|63| ) on (az — ni 2 ), and (120), (11S) and (|ll 



e 3 (a 2 + n 3 i) = e 2 (a 3 - rii 2 ) 



(56) 



There is still tetrad freedom remaining, but it is not obvious how one could best utilise it. The commu- 
tation functions which remain non-zero and the curvature variables they are coupled to in hi are &i, 2 , 
a Q , n 31 , n 12 , /j, and E n . This may be recognised as the Szekeres class I of dust solutions |23| (and its 
subcase, Ellis' dust spacetimes of LRS class II 0). The evolution equations and remaining constraint 
equations of these quantities arc also given there. 



B] ei(Hn) = 0: Immediately we see from the ei(f/n)-constraint (129) that ai Hn = 0, which provides 
us with the subdivision 

Bl] ai = or B2] H n = . 



Bl] oi = 0: We still have @ and @ holding, i.e., (a 2 - n 31 )if la = (a 3 + n u ) H n = 0. Now if 
Hn = 0, then we are dealing with B2]. So instead here we must have (a 2 — n^i) = (a3 + ni 2 ) = 0. We 
may use the tetrad freedom to set (03 — ni 2 ) = 0. From ( [75] ) and (115) we have (03 — ni 2 ) = without 
any further new constraints. And now, from (118), we have (03 — ni 2 ) = without any further new 
constraints. So this choice is allowed. Thus we may state that 



a.3 = n i2 = and (a 2 + n 31 ) = 2 a 2 . 

Now from ( |95|) combined with (|94|), ( |95| ) combined with (|96|), the commutator ( |6l| ) operating on //, and 
then using (pi|), we get that e a (/i) = e q(i^i i) = Al so, from ( |112| ) and ( |ll3| ) we have e 2 (9i) = e 2 (6 2 ) = 
^ e 2 (6) = e 2 ((Tn) = 0, while from (pT^ ) and ( JTl5| ) follows e 3 (6i) = e 3 (6 2 ) = ^> e 3 (6) = e 3 (o- u ) = 
0, which then, from the commutator (|6l|) acting on 9 and an, gives e a (9) = e a (<Jn) = 0. Moreover, it 
is clear from (|J) and (f|) that e a (Hn) = 0. The constraints (p6|), dn|) and ((ml) imply e a (n n ) = 0. 
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By applying the commutators (|5q), (59) and (|6C|) on the variables / S {ii, -En, 6, un, }, and 

), © 



utilising their respectiv e ev olution equations fl64|), (|79|), (|S8|), (|65|), ( |66| ) and (|76|), we can show that 
e a(f) — 0. Now from ( |12C| ) we have ei(et2) = 0, which allows us to use the remaining freedom to set 
e 3 (a 2 ) = 0. We can do this because 



e (e 3 (a 2 )) = - 2 6 2 e 3 (a 2 ) (horn fl6C|) on a 2 , (\[§ and (|U5|) ) , 
ei(e 3 (a 2 )) = (from © on a 2 and (|l20|) ) , 



and 

e 2 (e 3 (oa ) ) = 6 a 2 e 3 {a 2 ) (from (§on« 2 , (jl2|) and (p0|) ) 



Finally we set e 2 (a 2 ) — at an event x l — c\ This can be done because firstly, the eo-derivative of 
e 2 (a 2 ) is given by a multiple of e 2 (a 2 ) (from ( |59| ) operating on 02 and (|74|)); secondly, the ei-derivative 
of e 2 (a 2 ) vanishes (from ( |63[ ) operating on a 2 and ( |120[ )); thirdly, the e2-derivative of 62(02) is given 
by a multiple of 62(02) (from taking e2 of ( |122| )); lastly, the e 3 -derivative of 62(02) vanishes (from 
( |6l| ) operating on a 2 ). The remaining non-zero commutation functions and curvature variables they are 



coupled to in U are /i, 0i, 2 , a 2 , n\\, En and H\\. The remaining constraints (g_0Tp , (123) and (122)) 
are all algebraic relations. Thus the presently obtained consistent solutions constitute the subclass of 
spatially homogeneous cosmologies within the dust LRS class II of Ellis . 



B2] Hn = 0: From the H-equation ( p8[ ) it follows that nn an = 0. Now if nn — then we are dealing 
with case A] , and if in U an = 0, these are the FLRW dust solutions and are well known. 

2.4.2 Summary 

We conclude that the only non-trivial irrotational dust cosmologies that are PLRS according to Definition 



2.1 are known. These are the Szckcrcs dust cosmologies [|23j, which are strictly PLRS, and Ellis' dust 
spacctimes in LRS class II [pj. Moreover, this last study demonstrates that local rotational symmetry 
of an irrotational dust cosmology results if, in an open neighbourhood hi of an event of A4, all tensors 
algebraically determined from the Riemann curvature tensor and all its covariant derivatives up to only 
second order are rotationally symmetric, generalising a corresponding result given in 0. This property 
may be traced back to the fact that a rotationally symmetric covariant derivative Viei, i.e. 

T2ii = e 2 • Viei = = e 3 • Viei = r 3 n , 

corresponds by (|6l) and (J6$) to 



= (a 2 - n 3 i) = (a 3 + n 12 ) , 
eliminating the Szekeres dust models from the above consistent set and reducing to the LRS solutions. 



Note that these conditions also apply when p ^ ^> u\ ^ (cf. subsection 2.2 



3 Conclusion 

The only consistent solutions for the general class of perfect fluid cosmologies with all tensors rotationally 
symmetric about the same spatial axis are known solutions. All of these solutions are either LRS or they 
belong to the Szekeres class of dust cosmologies. Thus our work may be viewed as a form of classification 
scheme for inhomogenous cosmologies which incorporates the Szekeres dust solutions. This may be 
contrasted with the scheme developed by Szafron and Collins, where restrictions placed on submanifolds 
achieve this result ^2], A main result of the present work is that for all spatially inhomogeneous 



perfect fluid cosmologies that are PLRS according to Definition 2.1 we have in U 

H = . 

This generalises similar results obtained for the Szekeres dust cosmologies in JHj], and for perfect fluid 
cosmologies in LRS class II in jl3| . It has also been demonstrated that perfect fluid cosmologies are LRS 
if tensors and their covariant derivatives up to second order are rotationally symmetric, generalising a 
similar result obtained in l|7| . 
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There are no PLRS perfect fluid cosmologies which are rotating. Also, the cases which have vanishing 
fluid shear are fairly trivial: they are either not cosmologies at all in our understanding or they are of the 
simple FLRW kind. This last result requires the assumption of a barotropic equation of state p = p(n) in 
the irrotational accelerating perfect fluid case (although it is probably possible to relax this requirement). 

We have seen that the only strictly PLRS solutions are the Szekeres dust cosmologies. In later work 
we intend to relax the requirements of Definition 2.1 to see what partial symmetry features result. One 
possibility is to keep cr and E rotationally symmetric about the same spatial axis, but leave H unrestricted 
such that a uniquely defined tetrad {e Q } does exist. An example of a spatially inhomogeneous perfect 
fluid cosmology with coinciding eigenvalues of cr and non-zero E and H is provided by the diagonal 
Abelian G 2 solutions with radiation equation of state by Senovilla [Eol . 
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Appendix 

A Evolution and constraint equations 

Restrictions imposed: All 1 + 3 ONT equations we list in the following are written out for cosmologies 
with a perfect fluid matter source that are subject to the specific restrictions that 

• a fluid-comoving 1 + 3 ONT { e a } = { u, e a } is chosen which diagonalises cr, 

• cr is rotationally symmetric in the e 2 /e3~plane; meaning, 022 = 033 (= — \ en), 

• u>2 = = 0, that is, u) || e\, 

• 1x2 = U3 = 0; equivalently, u || ei. 

For brevity we use the variables 0i and 2 in some sets of equations while we use cru and O in other 
areas where the degeneracy in the eigenvalues of cr allows for some simplification. The relations 

011 = § (0i - 0a) and 6 = 6i + 2e 2 

yield the transformation from the former to the latter set of variables. 

The 1 + 3 ONT equations for perfect fluid cosmologies in fluid-comoving description without the 
specialisation stated can be accessed online at the URL O] . 



A.l Rotational tetrad freedom 

For coinciding eigenvalues of the fluid shear tensor cr such that 022 = 033 , its eigentetrad { e a } is defined 
only up to spatial rotations A = A(x l ) in the e2 / e3~plane. These may be represented by 



A" 



(1 

1 



\ 



\ 



cos ip sin tp 

- sin ip cos ip ) 



(57) 



with, in general, eo(v?) ^ 7^ e±((p). The members of the er-eigentetrad then transform according to 
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A. 2 Commutators 

In explicit form the commutators read 

[eo,ei](/) = uieo(/)-eiei(/)-n 3 e2(/) + n2e3(/) (58) 

[eo, e 2 ](/) = n 3 ei(/)-e 2 e2(/)-(wi+ni)e3(/) (59) 

[eo, e 3 ](/) = -n 2 ei(/) + (wi + ni)e2(/)-e 2 e3(/) (60) 

[e 2 , e 3 ](/) = 2wie (/) +nuei(/) - (a 3 -m 2 )e 2 (/) + (a 2 +n 3 i)e 3 (/) (61) 

[e 3 , ei](/) = (a 3 + m 2 )ei(/) + n 22 e 2 (/) - (ai -n 23 )e 3 (/) (62) 

[ei, e 2 ](/) = - (&2 - W3i)ei(/) + (ai +n 23 )e 2 (/) +n 33 e 3 (/) . (63) 

A. 3 Evolution equations 
A. 3.1 Energy density evolution 

Conservation of energy requires 

eo(/i) = -0(A* + p) ■ (64) 
A. 3. 2 Fluid shear and expansion evolution in terms of E 

floioi = -(eo + 9i)(ei) + (ei+ui)(tii) = Eu + ±(ji + 3p) (65) 
i?0202 = - (e + 2 ) (e 2 ) + w? - (oi + n 23 ) ui = £ 22 + ± (u + 3p) ( 66 ) 
i?0303 = -(e + e 2 )(6 2 )+a;f-(ai-n 2 3)wi = S 33 + i(/i + 3p) . (67) 

A. 3. 3 Fluid expansion evolution 

The expansion evolution is obtained by writing out Ro a oa- 

eo(0) - - | e 2 + ( ei + U! - 2 oi) («i) - § a\ x + 2 uo\ - § + 3p) . (68) 



A. 3. 4 Fluid vorticity and spatial commutation function evolution 



(023) 


e (a>i) = 


-2 9 2 wi - \ nn iii 


(69) 


(012) 


e (ai+n 23 ) = 


- 61 (ai + n 23 ) - (ei + iii) (6 2 ) + (e 2 - a 2 + 7131) (fi 3 ) 
+ (uji + fli) (n 22 - n 33 ) + 2 (»3 - n-12) 


(70) 


(013) 


e (ai-n 23 ) = 


- 61 (ai - n 23 ) - (ei + iii) (9 2 ) - (e 3 - a 3 - n 12 ) (0 2 ) 
- (wi + fii) (n 22 - n 33 ) - fi 3 (a 2 + n 3 i) 


(71) 


(012) 


e (a 2 -n 3 i) = 


- 6 2 (a 2 - n 3 i) - e 2 (©i) - (ei + «i - ai - n 23 )(fi 3 ) 
- (u>! + fii) (a 3 + m 2 ) - fi 2 (n 33 - nn) 


(72) 


(013) 


e (a 3 +ni 2 ) = 


- ©2 (ci3 + n 12 ) - e 3 (8i) + (ei + in - ai + n 23 ) (fi 2 ) 
+ (wi + fii) (a 2 - n 3 i) + fi 3 (nn - "-22) 


(73) 


(023) 


e (a 2 +n 3 i) = 


- (e 2 + a 2 + n 3 i) (6 2 ) + (e 3 - a 3 + n i2 ) (wi + fii) 
+ fi 3 (ai - n 23 ) + fi 2 (n 33 - nn) 


(74) 


(023) 


e (a 3 -ni 2 ) = 


- (e 3 + a 3 - n 12 ) (6 2 ) - (e 2 - a 2 - n 3 i) (wi + Qi) 
- fi 2 (ai + n 23 ) - fi 3 (nn - n 22 ) 


(75) 


(023) 


eo(nn) = 


- (i 9 - 2 (Tn) mi - (e 2 - 2 n 3 i) (fi 2 ) - (e 3 + 2 n 12 ) (fi 3 ) 


(76) 


(013) 


e {n 22 ) = 


- 61 n 22 - (ei + ui + 2 n 23 ) (wi + fii) - (e 3 - 2 n i2 ) (fi 3 ) 


(77) 


(012) 


eo(n 33 ) = 


- 6i n 33 - (ei + «i - 2 n 23 ) (wi + fii) - (e 2 + 2 n 31 ) (0 2 ) . 


(78) 
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A. 3. 5 Evolution equations for E and H 



[023] 23 



[031] 31 



[012] 12 



[031] 12 



[012] 31 



[012] 23 



[023] 12 



[031] 23 



[023] 31 



[023] 01 



[031] 02 



[012] 03 



[012] 02 



(e + 30 2 ){E 11 ) = (e 2 - a 2 - n 31 )(H 31 ) - (e 3 -a 3 + n 12 )(H 12 ) 

-±( f i + P )<j 11 + 2n 2 E 31 -2n 3 E 12 

-fniii/n + \{n 22 - n 33 )(H 22 - H 33 ) + 2n 23 H 23 

(e + ®){E 22 ) = -(ei + 2«i - en + n 23 )(H 23 ) + (e 3 - a 3 - ni 2 )(i?i 2 ) 

+ + p)trn + §eni£ 33 + 2f7 3 S 12 - (wi + 200^23 
-§n 22 ff 22 + |(n 33 - nn)(i? 33 - i?n) + 2n 31 H 31 

(e + 6)(-E 33 ) = (ei + 2iii - en - n 23 )(H 23 ) - (e 2 - a 2 + n 3 i)(fl" 3 i) 



(79) 



(80) 



+ + + 1^11^22 - 2f7 2 £ 3 i + (wi +20i)S 



23 



-|n 33 i? 33 + |(nn - n 22 )(Hu - H 22 ) + 2n 12 H 12 (81) 

(e + 6 + fcrnX^g) = -(ei + in - en + n 23 )(# 33 ) + (e 3 - 2a 3 - 2m 2 )(i? 3 i) 

+ i(M + pVi -ui(E 33 - E n ) - Cti(E 33 - E 22 ) - n 2 E 12 + n 3 E 31 

~2-(3n 22 + n 33 - nn)H 23 - (a 2 + n 31 )H 12 - (en - n 23 )H n + iiiH 22 (82) 

(e + 6 + \o xl )(E 23 ) = (ei + ui - ai - n 23 )(H 22 ) - (e 2 - 2a 2 + 2n 31 )(H 12 ) 

+ + wi(£; 22 - .En) - fti(£ 33 - e 22 ) + n 3 E 31 - n 2 E 12 

-|(3n 33 - nn + n 22 )H 23 + (a 3 - n 12 )H 31 + (en + rc 23 )i2ii - Mii? 33 (83) 
(eo + 6)(£?3i) = (ei + ui - 2oi - 2n 23 )(H 12 ) - (e 2 - a 2 + n 31 )(H n ) 
-n 2 (E n - E 33 ) - n 3 E 23 + (2wi + Oi)^i2 

-§(3rc 33 + nn - n 22 )H 31 - (a 3 + n 12 )H 23 - (a 2 - n 31 )H 22 (84) 
(e + 3e 2 )(£ 3 i) = (e 2 - a 2 - n 31 )(H 33 ) - (e 3 - 2a 3 + 2n 12 )(H 23 ) 
—Q 2 (Eii — E 33 ) + (ili — loi)Ei 2 — Q 3 E 23 

-|(3nn - n 22 + n 33 )H 31 + (ui + ai - n 23 )H 12 + (a 2 + n 31 )H 22 (85) 
(e + 6)(£i 2 ) = -(ei + tii - 2oi + 2n 23 )(£f 3 i) + (e 3 - a 3 - ni 2 ){Hn) 

-q 3 (£ 22 - En) + r* 2 £ 23 - (2wi + ni)s 3 i 

-i(3n 22 - n 33 + nii)i?i 2 + (a 2 - n 31 )H 23 + (a 3 + n 12 )H 33 (86) 
(e + 36 2 )(^i 2 ) = (e 2 - 2a 2 - 2n 31 )(H 23 ) - (e 3 - a 3 + n 12 )(H 22 ) 
—Q, 3 (E 22 — En) + (wi — Oi)^ 3 i + f2 2 -E 23 

-|(3nn + n 22 - n 33 )Hi 2 - (ui + ai + n 23 )H 3 i - (a 3 - n 12 )£f 33 (87) 
(e + 36 2 )(i?n) = -(e 2 - a 2 - n 3 i)(S 3 i) + (e 3 - a 3 + ni 2 )(£i 2 ) 

+ f nn^ii - i(n 22 - n 33 )(£ 22 - £ 33 ) - 2n 23 £ 23 + 2ft 2 ff 3 i - 2Q 3 H i2 (88) 
(e + 0)(H 22 ) = (ei + 2in - en + n 23 )(S 23 ) - (e 3 - a 3 - ni 2 )(£i 2 ) 

+ §e r ii J ff 3 3 + |n 22J B 22 - \{n 33 - n n )(E 33 - En) 

-2n 3 iE 3 i - (20i + Ui)H 23 + 2fl 3 Hi 2 (89) 
(e + ®)(H 33 ) = -(ei + 2iti - ai - n 23 )( J B2 3 ) + (e 2 - a 2 + n 3 i)(S 3 i) 
+ %a n H 22 + ^n 33 E 33 - \{nu - n 22 ){En - E 22 ) 

-2m 2 Ei 2 + (2fii + lui)H 23 - 2n 2 H 3 i (90) 
(e + 6 + §<7n)(# 23 ) = -(ei +u 1 -ai- n 23 )(E 22 ) + (e 2 - 2a 2 + 2n 31 ){Ei 2 ) 

+uJi{H 22 - Hn) - fti(H 33 - H 22 ) + Q 3 H 3 i - Q 2 Hi 2 - iei(n) 

+ 2-(3n 33 - nn + n 22 )E 23 - (a 3 - n 12 )E 3 i - (ai + n 23 )En + iiiE 33 (91) 
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[031] 03 



[012] 01 



[023] 03 



[031] 01 



[023] 02 



(e + 9 + f <jii)(# 23 ) = (ei + iii ~ ai + n 23 )(E 33 ) - (e 3 - 2a 3 - 2n 12 )(E 31 ) 

-ui(H 33 - H n ) - fii(£T33 - H22) - n 2 H 12 + fl 3 H 31 + ±ei(/i) 

+ i(3n 22 + n 33 - nn)E 23 + (a 2 + n 31 )E 12 + (ai - n 23 )En - iiiE 22 (92) 

(e + &)(H 31 ) = -(ei 2a x - 2n 23 )(E 12 ) + (e 2 - a 2 + n 31 )(E n ) 

-n 2 (H n - H 33 ) - n 3 H 23 + (2wi + fii)#i 2 + |eaO*) 

+ 5( 3 "33 + ™n - "- 22 )£ ; 3i + («3 + ^12)^23 + (a 2 - n 31 )E 2 2 (93) 
(e + 39 2 )(iJ 31 ) = -(e 2 - a 2 - ri 31 )(£ 33 ) + (e 3 - 2a 3 + 2n 12 )(£ 23 ) 
-n 2 (H u - H 33 ) - n 3 H 23 - (wi - Oi)#i 2 - ±e 2 (/x) 

+ i(3nn - n 22 + n 33 )-E 3 i - (tti + ai - n 23 )£i 2 - (a 2 + n 3 i)E 22 ( 94) 

(e + 6)(i? 12 ) = (e x +«i - 2ai + 2n 23 )(^ 31 ) - (e 3 - a 3 - n 12 )(£„) 
™Sl 3 (i? 22 - JT U ) + Sl 2 i/ 23 - (2w x + fixjffai - |e 3 (/x) 

+ i(3n 22 - n 33 + nu)£7i2 - (a 2 - n 3 i)£ 23 - (a 3 + rti 2 )£ 33 (95) 
(eo + 39 2 )(F 12 ) = -(e 2 - 2a 2 - 2n 3 i)(£ 2 3) + (e 3 - a 3 + n X2 )(E 22 ) 
-n 3 {H 2 2 - H u ) + n 2 H 23 + (ui - Sli)H 31 + 

+ |(3nn + n 22 - n 33 )E l2 + (ui + a x + n 23 )E 31 + (a 3 - n 12 )£; 33 . (96) 



A. 4 Constraint equations 

A. 4.1 Momentum conservation constraints 

We contract the second Bianchi identity (0) twice and find that for PLRS cosmologies according to 



Definition 2.1 



e i(p) + (t* + P)ui = 

62 (p) = 

es(p) = 



(97) 
(98) 
(99) 



A. 4. 2 Fluid vorticity, spatial commutation function and fluid shear constraints 



-R0303 


= (E 22 - E 33 ) = - 2 n 23 it! 






(100) 


(023) 


= E 23 = \ (n 22 - n 33 )ui 






(101) 


-Roio 2 


= E\i = e 2 (ui) +CL 2 u>i + |Sl 3 crn 






(102) 


-^0103 


= E u = e 3 (ui) +£l 3 ui - | n 2 crn 






(103) 


(012) 


0=| e 2 (?ii) - | (a 2 - n 3 i)ui + 2 u>i 






(104) 


(013) 


0=| e 3 (iii) - \ (a 3 + nu)ui + 3 wi 






(105) 


(123) 


= (ei — iii — 2 01) (wi) 






(106) 


-R0123 


= |rin (T11 + 2(tti + ai)o»i = -H n 






(107) 


-^0231 


= - (e x - ai + n 23 ) w x - | nn en - | (n 22 - 


n 33 )<Jn = 


- H 22 


(108) 


-^0312 


= - (ei - ai - n 23 ) cj x - § nn an + § (n 22 - 


n 33 )(7n = 


— H 33 


(109) 


-^0331 


= ei(6 2 ) + | (ai - n 23 )crn + 5 (n n + n 22 - 


n 33 )u>i = - 


- H 23 


(110) 


-^0212 


= -ei(6 2 ) - | (ax + n 23 ) a n + \ (n 22 - n 33 


-nn)u)x = 


- H 2 3 


(111) 


-Ron 2 


= e 2 (0i) - | (a 2 - n 3 i)a-n + (a 3 + ni 2 )a;i = 


= — H 31 




(112) 


-^0323 


= -e 2 (6 2 ) +e 3 (wi) = -H 31 






(113) 


-^0131 


= - e 3 (6i) + | (03 + ni 2 ) an + {a 2 - n 3 x)u>i 


= — H12 




(114) 


-Rq 2 23 


= e 3 (9 2 ) + e 2 (wi) = -H\ 2 






(115) 
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(123) 



R 



3112 



(123) 



R 



2312 — 



(l23) 



R 



2331 



R 



2323 — 



R 



3131 



i? 



1212 



ei(nn) + e 2 (a 3 + n 12 ) - e 3 (a 2 - n 3i ) - 2a x n u - 2a 2 n i2 - 2a 3 n 3i 
= -26iwi 

- e 3 (a 2 - n 3 i) + | ei(rin + n 22 - n 33 ) - nu (01 - n 23 ) 

- ai ("22 - "-33) - "-23 ("-22 + n 3 3) - 2 n 3 i (a 3 + n i2 ) + 61 wi 

- B 23 

e 2 (« 22 ) - e!(a 3 - n 12 ) + e 3 (ai + n 23 ) - 2a x n 12 - 2a 2 n 22 - 2a 3 n 23 
= — 2 3 wi 

- e 3 (ai + n 23 ) - | e 2 (nn + n 22 - n 33 ) + nu (a 2 - n 3 i) 
+ (a 2 + n 31 ) (n 22 - n 33 ) +2n 23 (a 3 - n 12 ) - 2fi 3 wi 

- -B31 

e 3 (n 3 3) - e 2 (ai - n 23 ) + ei(a 2 + n 3 i) - 2ai n 3i - 2a 2 n 23 - 2a 3 n 33 
= 2 f2 2 wi 

- e 2 (ai - n 23 ) - \ e 3 (n 22 - n 33 - nu) - nu (a 3 + n X2 ) 
+ (a 3 - m 2 ) (n 22 - "33) ~ 2n 23 (a 2 + n 3i ) -2fi 2 wi 

- i? 12 

e 2 (a 2 + n 3 i) + e 3 (a 3 - n i2 ) + 6 2 + cj 2 

- (ai + n 23 ) (ai - n 23 ) - (a 2 + n 3 i) 2 - (a 3 - n i2 ) 2 

- I "11 + 5 "11 ("22 + "33) + 2 (n 2 2 - n 33 ) 2 - 2 fii uj t 



e 3 (fl3 + "12) + ei(ai - n 23 ) + 6 2 61 

- (a 2 + n 3 i) (a 2 - n 3i ) - (a 3 + n i2 ) 2 - (ai - n 23 ) 2 

+ i n ii + i ("-22 - "-33) (2 nu - 3 n 22 - n 33 ) 



= -J5 



22 



ei(ai + n 23 ) + e 2 (a 2 - n 3i ) + 6 2 61 

- (a 3 + m 2 ) (a 3 - n X2 ) - (ai + n 23 ) 2 - (a 2 - n 31 ) 2 



1 T7 2 

4 n n 

Br- 1 1 



"T 1 - \ (n 22 - n 33 ) (2 nu - 3 n 33 - n 22 ) 
:;.-) t- 3 M • 



(116) 

(117) 
(118) 

(119) 
(120) 

(121) 
(122) 
(123) 
(124) 



A. 4. 3 Generalised Friedmann equation 



This is obtained by writing out R al3 a /3 (which corresponds to summing ( 122 ), (|l23| ) and (124)) 



4ei(ai) + 4e 2 (a 2 ) +4e 3 (a 3 ) - 6a 2 - 6aj -6a 2 . 



2 n 23 



2 "li 



2 P)2 _ 3 2 



2n\ 2 



in 2 



11 



I (n 22 - n 33 ) 2 + nu (n 22 + n 33 ) 



2w 2 - 4Qi Wi = 2/i 



(125) 



A. 4. 4 Constraints on E and H 



[123] 23 



[123] 31 



[123] 12 



[123] 01 



(ei - 3ai)(E n ) + (e 2 - 3a 2 + n 31 )(E 12 ) + (e 3 - 3a 3 - ni 2 )(E 31 ) 
iei(/i) + (n 22 - n 33 )B 23 - (E 22 - E 33 )n 23 + 3wiH n 



(e 2 - 3a 2 )(i? 22 ) + (e 3 - 3a 3 + n 12 )(E 23 ) + (ei - 3ai - 
|e 2 (/x) + (n 33 - nu)E 3 i - (E 33 - E n )n 3 i + 3u>iH 12 
(e 3 - 3a 3 )(^ 33 ) + (ei - 3a x + n 23 )(£ 3 i) + (e 2 - 3a 2 - 
f e 3 (/z) + (nu - n 22 )£a 2 - (En ~ E 22 )n 12 + 3u>iH 31 
(ei - 3ai)(ffn) + (e 2 - 3a 2 + n 3 i)(ifi 2 ) + (e 3 - 3a 3 - n 12 )(H 31 ) 



n 23 )(Ei 2 ) 
- 2 &iiH 3 i 
n 31 )(E 23 ) 



\<J\\H\ 2 



(126) 
(127) 
(128) 
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[123] 02 



[123] 03 



- (/i + p)u>i + (n 22 - n 33 )H 23 - (H 2 2 - H 33 )n 23 - 3wiEn (129) 

(e 2 - 3a 2 )(H 22 ) + (e 3 - 3a 3 + rai2)(i?23) + (ei - 3ai - 7123) (ifo) = 

(n 33 - nu)H 31 - (H 33 - Hn)n 31 - "iujiE 12 + |(Tii-B 3 i (130) 

(e 3 - 3a 3 )(if 33 ) + (ei - 3ai + n 23 )(H 31 ) + (e 2 - 3a 2 - n 3 i)(i?23) = 

(nn - n 22 )H 12 - (Hu - H 22 )n 12 - ZujiE 31 - \o\\E X2 . (131) 
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